Normal and lateral critical Casimir forces between colloids and patterned substrates 
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We study the normal and lateral effective critical Casimir forces acting on a spherical colloid immersed in 
a critical binary solvent and close to a chemically structured substrate with alternating adsorption preference. 
We calculate the universal scaling function for the corresponding potential and compare our results with recent 
experimental data [Soyka R, Zvyagolskaya O., Hertlein C, Helden L., and Bechinger C, Phys. Rev. Lett., 101, 
208301 (2008)]. The experimental potentials are properly captured by our predictions only by accounting for 
geometrical details of the substrate pattern for which, according to our theory, critical Casimir forces turn out to 
be a sensitive probe. 

PACS numbers: 

05.70.Jk — Critical point phenomena 
82.70.Dd — Colloids 

68.35.Rh — Phase transitions and critical phenomena 



Introduction. — The confinement of a fluctuating 
medium generates effective forces acting on its boundaries. 
A particularly interesting realization of this general princi- 
ple is provided by the confinement of concentration fluctua- 
tions of a binary liquid mixture upon approaching a critical 
demixing point at temperature T = T c in its bulk phase dia- 
gram [1]. Generically, the confining surfaces preferentially 
adsorb one of the two components of the binary liquid. This 
amounts to the presence of effective, symmetry-breaking sur- 
face fields favoring either positive [(+)] or negative [(— )] val- 
ues of the scalar order parameter which is the difference 
between the local concentrations of the two species. The ex- 
tension of the spatial region in the direction normal to the 
surfaces, within which the local structural properties of the 
fluid deviate from the bulk ones, is given by the bulk corre- 
lation length 4 , which diverges upon approaching the critical 
point as § (t -> 0) = §f \t | ~ v . Here t = (T — T c ) /T c is the re- 
duced temperature 1 , v ~ 0.63 is a standard bulk critical expo- 
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Hi t = {T c - T)/T c . 



nent, and §q are nonuniversal amplitudes for t ^ 0. In a film 
of thickness L, the ensuing long-ranged critical fluctuations 
of the order parameter lead to a critical Casimir force acting 
on the confining walls (J], which is described by a universal 
scaling function of Lj% and depends on the type of boundary 
conditions (BC) at the walls (see, e.g., Refs. Jj, |5l la |7|] and 
references therein). This is the thermodynamic analogue of 
the quantum-electrodynamic Casimir effect originating from 
the confinement of vacuum fluctuations |8l|9|]. Depending on 
the relative adsorption preferences of the boundaries, the crit- 
ical Casimir force is either attractive for identical symmetry- 
breaking BC (±,±) or repulsive for opposite BC (±,=f). 
(Symmetry-preserving BC can lead to attractive and repulsive 
critical Casimir forces, too fl4 ll011 .) Besides various indirect 
experimental evidences in thin films [ 1 1], the critical Casimir 
force has been measured also directly at the sub-micrometer 
scale for a colloid immersed in a near-critical binary liquid 
mixture close to a homogeneous substrate [2]. Recent Monte 
Carlo simulations are in quantitative agreement with all avail- 
able experimental data Jjj [12J, ll3|] . 

Colloids can be used not only as model systems in soft mat- 
ter physics but also in applications on the nano- and microme- 
ter scale which take advantage of their interaction with chem- 
ically structured solid surfaces. Such systems can be useful 
in integrated nano-devices provided that one is able to exert 
active control over these interactions. Critical Casimir forces 
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provide such a tool, because their strength and direction can 
be tuned via minute temperature changes and surface treat- 
ments of the substrate. Recently, the critical Casimir potential 
of a colloid close to such a chemically patterned substrate has 
been measured [3], providing evidence for the occurrence of 
lateral critical Casimir forces. In Ref. Jst], a dilute suspension 
of charged spherical colloids, imposing (— ) BC to the order 
parameter of the near-critical solvent, is exposed to a chemi- 
cally patterned substrate, the surface of which consists of al- 
ternating stripes, which impose (— ) and (+) BC. The equilib- 
rium spatial distribution of colloids was measured via digital 
video microscopy and from it one can define an effective po- 
tential for a single colloid, which varies laterally due to criti- 
cal Casimir forces. In view of potential applications, surfaces 
might be designed as to provide temperature-controlled lat- 
erally confining potentials for single colloids, offering novel 
means of self-assembly processes yfl. 

In order to reap the full benefits of this wide range of pos- 
sibilities, a thorough theoretical understanding of the underly- 
ing physics is essential. In this context, lateral critical Casimir 
forces are theoretically known to occur in the film geometry 
with chemically [14] or geometrically lfl5ll structured sub- 
strates. However, in the experimentally relevant realizations 
of such a film geometry, based on wetting phenomena |4j,|9|] at 
least one of the two confining surfaces (i.e., the liquid-vapor 
interface) is laterally homogeneous and therefore the lateral 
critical Casimir force vanishes. Instead, for the geometry of 
a colloid facing a substrate such a lateral force is expected to 
arise even if only the substrate is laterally inhomogeneous, as 
the experimental findings in Ref. ]3|] demonstrate. Theoret- 
ical studies for the geometry of a colloid facing a wall have 
been limited to laterally homogeneous surfaces lfl6l [nil . In 
order to overcome these limitations and to interpret, inter alia, 
the experimental data of Ref. yfl, we have studied the critical 
Casimir force acting on a sphere of mesoscopic radius R at 
surface-to-surface distance D from a substrate with laterally 
alternating adsorption preferences for the two components of 
a confined binary liquid mixture. We provide quantitative pre- 
dictions for the universal features of this effective force (i.e., 
in excess to regular, nonuniversal background contributions) 
by using the so-called Derjaguin approximation (DA) together 
with the knowledge of the scaling functions determined by 
Monte Carlo simulations in spatial dimension d = 3 for the 
laterally homogeneous film geometry IU2I1 . In order to esti- 
mate the accuracy of this Derjaguin approximation, we per- 
form a full numerical analysis of the appropriate mean-field 
theory (MFT) without further approximation and compare it 
with the corresponding Derjaguin approximation based on an- 
alytic results in d = 4. 

The demixing point is approached from the mixed phase 
by varying the temperature T towards T c at fixed pressure 
and critical composition of the binary fluid. Within the field- 
theoretical renormalization group approach the fixed-point 
Hamiltonian for a binary mixture is given by 11811 

^]=^d d r|I(V^) 2 + ^ 2 + ^ 4 |, (1) 

where the integration runs over the volume V accessible to the 



fluid described by the af-dimensional position vector r, T«(, 
and u > is a coupling constant. In the strong critical ad- 
sorption limit lfl9tl the surface contributions to the Hamilto- 
nian turn into BC corresponding to infinite surface fields so 
that 0| surfa e = ±°°. Thus, J^[<j)] is supplemented by the BC 
= — °o [{—)] at the surface of the colloid and by (j) = — °° 
(+oo) on that part of the substrate with the same (opposite) 
adsorption preference. Using a 3d finite element method, we 
have numerically minimized J^[(p — and have ob- 

tained the mean-field order parameter profile m. This MFT 
solution allows one to infer the universal scaling function of 
the critical Casimir force at the upper critical dimension d = 4 
up to an overall prefactor °< u (and up to logarithmic cor- 
rections). The normal and lateral critical Casimir forces are 
calculated using the stress tensor l20ll and the associated po- 
tential is determined by integration. 



Homogeneous substrate. — We first consider a three- 
dimensional sphere with BC (b) facing a homogeneous sub- 
strate with BC (a), denoting this combination as (a,b). 
The critical Casimir potential &( a ,b) takes the scaling form 2 
® M (D,R,T)=k B TRD 2 - d $ (atb) l®,A), where A = D/R, 
© = ±D/£± for t ^ 0, and 1?^) is a universal scaling 
function yl [13]. If D <C R, i^L&ican be expressed - 
via the Derjaguin approximation [T3, I2TI1 - in terms of 
the critical Casimir force per unit area ff a ^(L,T ^ T c ) = 
kgT L~ d k( a b }(±L/ 1;±) between two planar walls at distance 
L with (a,b) BC. For d — 3,4 the Derjaguin approximation 

gives U\m 

(fli6) (0,A^O) =2k £ds(s 1 - d -s- d )k M (s&). (2) 

Comparing in d = 4 the numerically calculated scaling 
functions to the approximate ones obtained via the 

Derjaguin approximation (on the basis of £(±.-) given in 
Ref. D20D), we infer that the former are reasonably well de- 
scribed by the latter for A < 1/3. Assuming a smooth depen- 
dence on d, we expect this to hold for d = 3 as well. 



Chemically patterned substrate. — The essential build- 
ing block of a substrate with laterally varying adsorption pref- 
erence is a chemical step, i.e., two semi-infinite planes with 
(ag) BC for x ^ joined together, where x is one of the lat- 
eral position coordinates; with respect to the latter the center 
of the spherical colloid with (b) BC is located at x = X. The 
resulting critical Casimir potential <t> for d = 3,4 can be cast 
in the form 

<S>(X,D,R,T) =k B TRD 2 - d #(S,@,A), (3) 



2 In d = 4 the 3d sphere is a hypercylinder and &( a ,b) is the potential per 
length in the fourth direction, along which the physical properties are in- 
variant. 



3 



where E = X j\jRD is the scaling variable associated with the 
lateral position of the colloid. For E — > ±°°, i.e., sufficiently 
far from the step, along the planar substrate one recovers the 
homogeneous cases, i.e., #(E — > ±°°,0,A) = #( a> m(0,A). 
The normal and lateral critical Casimir forces acting on the 
colloid are given by — dp<P and — dx$>, respectively. It is con- 
venient to introduce the scaling function co defined via 

0(3,0, A) = *("<•»)+*'''>•») + gfa^Q^a fl) (B,Q,A) 1 (4) 

where fy a> ,b) = ^(a>,fc)(®,A) re f er to homogeneous 
substrates and depend only on © and A, while 
©(E — > ±°°,0,A) = Note that the common prefac- 

tor °c k , which within MFT is left undetermined by the 
analytical and numerical determination of t?( fl> &)(©, A) 
and #(S,0,A), does not affect the MFT prediction for 
<o(E,0,A). 

For small distances D <^R, i.e., A — > 0, the surface of the 
colloid facing the substrate can be considered to be made 
up of successive and consecutive circular rings, parallel to 
the substrate, which have an infinitesimal area dS(l) and a 
radius r(l) which increases with increasing the normal dis- 
tance I from the identical corresponding circular ring obtained 
by normally projecting the ring on the surface of the sphere 
onto the one of the substrate. Assuming additivity the criti- 
cal Casimir force F acting on the colloid is the result of the 
sum of the forces dF(l, . . .) acting on each single pair of such 
corresponding rings with separation /. Neglecting also edge 
effects, dF(/,...) can be expressed in terms of the critical 
Casimir force per unit area fi a> ,b) (hT) acting on infinite and 
homogeneous parallel plates with (a^,b) BC separated by the 
same distance / as the rings: dF(l, ■■■) — fi a> b){^T)dS > + 
f {a<ib) (l,T)dS<. Here dS% (with dS(l) = dS> +dS<) indi- 
cates the surface area of that portion of the ring on the sub- 
strate for which x ^ 0, i.e., corresponding to (a^) BC. This 
area dS^ depends, inter alia, on / and on the position of the 
colloid. For the potential <t> (see Eq. (01) associated with F 
we eventually find in d = 3 , 4 

0(S % 0, 0, A - 0) = # M (0, A -> 0) ± IE 4 

t s arccos (s~ 1 / 2 ) — Ws — 1 , . 1 , (5) 

x/ds i '—V A/t(0(l+E 2 s/2)), 

/ (l+Zh/2) d V V " 

where Ak(&) — k/ a< m(0) — &( a> ,*>)(©) is me difference be- 
tween the scaling functions of the critical Casimir forces 
acting on two planar, homogeneous walls with (a <1 b) and 
(a > ,b) BC, respectively. Equation (0 yields 0)(E = 0,0, A — > 
0) = (expected from Eq. and additivity) and 

O)(E,© = 0,A^0) =S(l-rf-E 2 )(E 2 + 2) _3/2 , (6) 

independent of k( a> b y Note that Eqs. d3j-(|6]l are valid beyond 
the BC we consider in the following, i.e., S {+,— } and 
b = — . For (=p, — ) BC the critical Casimir force /(=p -) (D, T) 
between two planar walls decays <^ exp(— 0) as a function 
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FIG. 1: (Color online) Scaling function CO [Eq. ©] within Derjaguin 
approximation (i.e., A — > 0) for the critical Casimir potential of a 
sphere (— ) facing a chemical step (+[—). The dashed and dotted 
lines refer to d = 3 and d = 4 (MFT), respectively. At bulk criticality 
(0 = 0), CO is given by Eq. ((6), whereas for ^ it is calculated 
on the basis of the scaling functions for the film geometry [see the 
main text]. For > 3, CO becomes practically independent of d and 
coincides with the expression for 2> 1 [Eq. l(7J, solid lines, barely 
distinguishable from the corresponding dashed and dotted ones]. 

of for 0> 1 12(1 |22l] 3 , which leads to a d-independent 
result for co: 

<o(S,0> 1, A -+0) =-erf( v /©72E) ) (7) 

where erf is the error function. 

In Fig. Q]we compare the behavior of CO calculated within 
the Derjaguin approximation in d = 4 and 3; the required scal- 
ing functions kt±_\ for the film geometry [see Eq. (0] are 
obtained analytically in d = 4 within MFT [20] and in d — 3 
from Monte Carlo simulations B12I1 . The systematic uncer- 
tainty of the latter does not affect significantly (at most by 
3%) the estimate of co shown in Fig.Q] For 0^0 the critical 
Casimir potential for d = 3 as a function of the lateral coordi- 
nate varies more smoothly than the corresponding MFT func- 
tion [see Eq. ©]. However, for > 3 the scaling functions co 
for d = 3 and d = 4 (MFT) practically coincide with Eq. (j7]) 
valid for 0^1. Figuref2]compares the scaling function co for 
d = 4 (MFT) calculated in the limit A — > (Derjaguin approxi- 
mation) with the full one determined numerically for A = 1/3. 
The former provides a very good approximation of the latter, 
especially for large 0. Since a viable strategy for carrying out 
Monte Carlo simulations for the pertinent lattice Ising model 
in the sphere-plate geometry is currently not available, this 
study of the DA in d = 4 is an important piece of information 
because at present this is the only way to quantitatively assess 



3 This purely exponential decay is proven analytically only in d = 2 and 4 
but it is expected to hold for all spatial dimensions. 
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FIG. 2: (Color online) Scaling function (0 [Eq. @J] in = 4 (MFT) 
for the critical Casimir potential of a sphere (— ) facing a chemical 
step ( + 1 — ) . The numerical results for A = 1 /3 are indicated by sym- 
bols, the size of which refers to the numerical error. The lines show 
the corresponding results for A — > 0, calculated within the Derjaguin 
approximation, which provides a good approximation to the numeri- 
cal data for > 3. 



the range of parameters within which the DA provides an ac- 
curate approximation of the actual scaling function. On the 
basis of these results in d = 4 and assuming a smooth depen- 
dence on d, we expect that also in d — 3 the DA captures the 
actual behavior of the scaling function for small A and large 
© — two conditions which are met by the experimental data 
discussed below. 

A consequence of the underlying assumption of additiv- 
ity is that the critical Casimir potential due to an arbitrarily 
striped substrate can be calculated within the Derjaguin ap- 
proximation on the basis of Eq. (0 as the appropriate super- 
position of consecutive, shifted steps. From the comparison of 
the Derjaguin approximation for this potential in d = 4 with 
the corresponding MFT results we find that the former de- 
scribes quite well the actual behavior predicted by the latter 
even for substrates with a fine pattern, at least as long as the 
ratios between the characteristic lengths of the pattern and the 
geometrical average \/RD of the radius R of the colloid and its 
distance D from the substrate are not too small. In particular 
for a chemical stripe of finite width L > 0.2\/RD the use of 
the Derjaguin approximation for A < 1/3 is justified. More- 
over, for Lj VRD > 3 and ® > 3 or for L/y/RD > 2 and © > 8 
the stripe is almost equally well described by two independent 
and subsequent chemical steps. Similarly for periodically pat- 
terned substrates of alternating stripes with total period P the 
DA describes the actual data accurately if, in addition to the 
fulfillment of the previous conditions for each single stripe in 
the pattern, one has P > 3\/RD. We expect that this property 
carries over to the case of d = 3. 



Comparison with the experiment by Soyka et al. — In 

Ref. [3J the substrate surface confining the colloidal suspen- 
sion consists of stripes along the y-axis which impose (— ) 
and (+) BC alternating along the x-axis, and have a width 
L_ = 2.6/im and L + = 5.2 fim, resulting in a periodicity 
P = L_+L+ =7. 8. ^m. 

The total potential 4> tot (in units of IcbT with T ~ T c ~ 
307 K) of the forces acting on each colloid (R = 1.2/xm, 
(— )BC) is the sum of the electrostatic, gravitational, and crit- 
ical Casimir contributions: 

*tot(jc,£>) =<S>d(x,D) + <t> g (D)+<t> c (x,D). (8) 

In general the colloidal particle interacts with the substrate 
and the surrounding medium also via van der Waals forces[2, 
l23ll . the potential of which should be added to the rhs of 
Eq. ([8]l. 4 However, it has been shown that for the particu- 
lar choice of materials and experimental conditions used in 
Ref. JH, the corresponding contribution is negligible com- 
pared with the others [IHO]. For the equilibrium number den- 
sity of colloids one has p(x,y,D) °< exp{— <3? tot (x,D)}, where 
x and y are the coordinates of the projection of the center of the 
colloid onto the substrate surface. Since the distance D is not 
resolved in the setup of Ref. yl, only the projected number 
density J^dD p(x,y,D) is experimentally accessible. Due to 
the translational invariance of the pattern along the y-direction 
(of length I 3> R,L±) one can in addition project the experi- 
mental data onto the x-axis, t~ 1 f dy J$dD p(x,y,D) =: p(x), 
and define, up to a constant, an effective potential V such that 
p(x) exp{— V(x)} J3|]. Our analysis shows that for the ex- 
perimental conditions used in Ref. [13] the effects of consec- 
utive chemical steps do not interfere and therefore one can 
focus on a single step located at x = 0, with (=p) BC for 
x ^ 0. At this stage we assume that each chemical step ex- 
perimentally realized on the substrate is effectively sharp and 
straight at the micrometer scale of the problem. For such 
a step, 4> e i(x,D) interpolates between the screened electro- 
static potentials Oei.^ (D) of a colloid facing a homogeneous 
substrate, which are recovered far from the step for x ^ 
and which are well approximated by exp{— (D — D^)/X} 
where A ~ 12nm is the screening length of the mixture and 
~0.1-0.2/imj2[]. Within the Derjaguin approximation 

4> el (x,D) =4> el ,+ (D)0(-x/A)+<I> eL _(D)0(x/A), (9) 

where 0(h) = [1 +erf(u)]/2 and A = VZRX ~ 0.17 jum. 
The gravitational potential in units of k%T at T ~ 307K 



Upon approaching the critical point, the temperature-dependent dielectric 
permittivity e(T) of the binary liquid mixture exhibits a weak cusplike 
singularity e(T) — e(T c ) oe |f|'~ a , where a ~ 0.11 is the critical exponent 
of the specific heat for the three-dimensional Ising universality class 1 24]. 
This weak variation of e(T) might affect the strength of the van der Waals 
forces, as well as the range of the electrostatic interaction 3>,.i 1 2]. However, 
in the near-critical mixture of water and lutidine, the permittivity e(T) turns 
out to vary less than 1% for \T C — T\ < IK I2al , and therefore the van der 
Waals forces as well as the electrostatic interaction are expected to be not 
affected significantly. 
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is %(D) = GD with G = (p PS - p v/L )(47zR 3 /3)g/{k B T) 



for a gravitational acceleration g 



9.8 m/s^ 



1.12 Gum)- 1 

where ppg ~ 1.055 g/cm J and Pwl — 0.988 g/crrr' [26] are the 
mass densities of the polystyrene colloid and the solvent, re- 
spectively. 4>c is calculated on the basis of Eqs. (f3]) and (|5]l 
with #(±,-) given by the Derjaguin approximation for ho- 
mogeneous substrates and k(±~) obtained from Monte Carlo 
simulations [12]. Far from the critical point <t>c is negligi- 
ble compared to <t> e i,g and the average height (D) of the col- 
loid above the substrate is ~ + G~ l ~ 1 fim, with typical 
fluctuations ~ G ~ 0.8/xm. (These values are significantly 
larger than those reported in Ref. [3], which were extrapolated 
from the different experimental conditions used in Ref. J2l .) 
Upon approaching the critical point, the repulsive critical 
Casimir force for x < pushes the colloid slightly away from 
the substrate, causing only a minor increase of (D) . For x > 0, 
instead, as soon as the correlation length 4 exceeds a certain 
Dq -dependent value, the colloid abruptly localizes in the po- 
tential well due to the interplay between the attractive critical 
Casimir force and the electrostatic repulsion, with (D) ~ Dq 
and fluctuations of few tens of nm. As a result, at distances 
D ~ (D), the contributions of 4> el to <£> tot (x < 0,D) and of 4> g 
to <i>tot{x > 0,D) are negligible and the behavior of the colloid 
depends sensitively only on the actual value of . Accord- 
ingly, we fix Dq =0.1 jxm, X = 12nm and calculate V(x) and 
AV := V (-L+/2) - V(L-/2) on the basis of <P tot (x,D) for the 
single step. V and AV depend on Dq and, via the bulk correla- 
tion length 4 ~ | + [(Ar c + AT) /307 K] -°- 63 , on the estimated 
distance AT = T c — T from the critical point (actually lo- 
cated at T = 307K + AJ C ). The values D ~ 0.136jUm, 4 + - 
0.42 nm, and AT C ~ 27 mK yield a very good fit to the experi- 
mental data for AV (Fig. 3 of Ref. Whereas AT C is within 
the experimental accuracy, 4^ is significantly larger than pre- 
vious estimates 4^ = (0.2 ± 0.02) nm ||2|], suggesting that for 
the rather small corresponding values of 4 — 20-36 nm cor- 
rections to the leading scaling behavior might still be relevant. 
For x < the contribution of 4>c to 4> tot is significant only 
for distances from the substrate D < 4 , which corresponds to 
A = D/R < 0.03, whereas for x > the typical distance is 
D~D and therefore A ~ 0.12 with © = Dj% > 4. In addi- 
tion, the relevant geometrical parameters are P/^/RD ~ 7-20 
and L-j\jRD ~ 2-7, such that, based on our theoretical anal- 
ysis of the range of validity of the DA, we expect the latter 
to be accurate for the potentials <J> e i and <t>c- Moreover, the 
values of ©, P/\/RD, and L-j\fRD are such that the result- 
ing critical Casimir potential is adequately described by the 
superposition of a sequence of single chemical steps. 

In Fig. [5] we compare the experimental data with the re- 
sulting theoretical predictions for 8V(x) :— V(x) — V(—L + /2) 
(solid lines) across the (— ) stripe (grey). As anticipated, the 
effects of the chemical steps at x — and L do not inter- 
fere in the actual range of parameters. The theoretical curves 
- in qualitative disagreement with the experimental data - 
display a sharp transition between the plateau values and 
—AV. In order to test the robustness of this distinctive fea- 
ture we varied X within the plausible range 8-1 8 nm, con- 
sidered polydispersity (R = 1.2-1. 8 jUm), and allowed for a 




FIG. 3: (Color online) Lateral variation of the effective potential 8V 
[see the main text] of a colloidal particle (— ) facing a chemically 
patterned substrate and immersed in a binary liquid mixture at critical 
concentration for various temperatures T c — AT. Symbols indicate 
the experimental data of Ref. |3], whereas the solid and the dashed 
lines are theoretical predictions for the same values of parameters for 
an ideal and non-ideal stripe pattern, respectively, based on Eq. ((5J. 



possible inhomogeneous buoyancy <P g (D,x) induced by lat- 
erally varying fluid layers adsorbed on the colloid and on 
the substrate. Moreover, due to their fabrication process, the 
(— ) stripes might have a rather weak preferential adsorption 
as compared to the (+) stripes, which we tried to capture 
by reducing the amplitude of -) by up to 70%. Conse- 
quently, the values of Dq ~ 0.10-0. 14jtim, 4 + — 0. 3-0.4 nm, 
and AT C ~ —0.1-0.1 K, which yield the best agreement with 
the experimental data for AV, are affected by these changes 5 
but the sharpness of the variation is not. In addition, our anal- 
ysis shows that the effect of the periodic spatial arrangement 
of the stripes - which would smoothen the potential for rela- 
tively small periodicities P - is negligible for the experimen- 
tal conditions used in Ref. |Ht]. On the same basis, we expect 
the DA to be sufficiently accurate and that the non-linearities 
inherent in the critical Casimir interaction, which actually in- 
validate the assumption of additivity of the forces and might 
cause a smoothening, do not affect significantly the potential 
shown in Fig. [3] 

However, if the actual position x = x s (y) of each chemi- 
cal step varies smoothly along the y-axis on the length scale 
I ~ 40 fim explored by the colloid during the measurement, 
the projection of the number density p onto the x-axis results 
in a smoothed distribution p(x). Even though there is no di- 
rect measurement of such a variation of the position of the 
step, it is reasonable to assume that it occurs on the length 



' As the amplitude of v is reduced, the fitted value of moves closer 

to Iq 1 ", suggesting that indeed the preferential adsorption of the (— ) stripes 
might be effectively rather weak. 
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scale I due to the fabrication process (focused ion beam act- 
ing on glass) and due to the projection of the data onto one 
dimension without independent knowledge of the alignment 
of the chemical stripes. In order to estimate the consequences, 
we have assumed the total potential of the forces to be given 
by <S>tot(x — x s (y),D) with x s (y) characterized by a Gaussian 
distribution p(x s ) with zero average and standard deviation 
Ax = 0.5 fim. (Ax may contain a contribution from a smooth 
intrinsic chemical gradient.) Thus, for the number density, the 
projection IT 1 § dy turns into f^™dx s p(x s ). The resulting 8V, 
which basically correspond to a convolution of the original al- 
most square-well-like potentials with a Gaussian of width Ax, 
are shown as dashed curves in Fig. [3] The agreement with the 
experimental data is significantly improved. In view of the 
pronounced modifications of the resulting potential our anal- 
ysis demonstrates that critical Casimir forces respond sensi- 



tively to geometrical details of the chemical pattern, which 
could not be checked independently in Ref . y|] . 

We conclude that our theoretical analysis provides a quanti- 
tative understanding of lateral critical Casimir forces and new 
insights into the corresponding experiments. The reliability of 
the theoretical description allows one to use critical Casimir 
forces for probing microscopic features of the system which 
are difficult to access otherwise. 
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